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ABSTRACT 



The design of a supersonic missile seeker lens system 
using an aerodynamically efficient lens, a scanning mirror 
and detector was accomplished through the use of ray- 
tracing routines. The design of the first element lens 
uses a gradient in refractive index to overcome the severe 
handicap that the pointed shape imposes. This problem has 
been previously solved for positive orientations of the 
center of symmetry of the Gradient Refractive Index (GRIN) 
material with respect to the lens. In this study the 
negative orientations are resolved and the lens studies for 
a range of orientations and strengths of gradient. The 
best lens is then used in conjunction with the mirror and 
detector in a ray- tracing scheme. Such a system was found 
to be feasible and worthy of further study. 
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I. 



INTRODUCTION 



A. OUTLINE OF LENS DESIGN TO DATE 

The design of a 'pointed' aerodynamically efficient 
missile seeker lens has been the topic of a number of 
studies. Frazier [Ref. 1], and Terrell [Ref. 2 ], first 
studied the idea of using a pointed lens as a means of 
obtaining some degree of aerodynamic efficiency. They also 
introduced the idea of using Gradient Refractive Index 
(GRIN) material in order to overcome the difficulties in 
imaging imposed by the odd shape. Amichai [Ref. 3], and 
Carr [Ref. 4], studied the GRIN concept in more detail, 
Amichai by defining the outer surface as a cone and solving 
for the inner surface and Carr by defining the inner 
surface as a cone and solving for the outer. The main 
reference for this author's study was drawn from Carr's 
work and a review of that thesis would be helpful if the 
reader seeks a more detailed background from that given in 
the next few pages. 

B. GRADIENT REFRACTIVE INDEX (GRIN) 

Although the theory for GRIN has been developed and 
known for some time it was not until fairly recently that 
materials have become available that allow the theory to be 
put into practice. A gradient in index will cause light 
rays to bend in an otherwise homogeneous material. Thus, 
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the lens designer now has a new dimension of flexibility 
which will allow him to solve problems that could not be 
formerly solved. A proper selection of the manner in which 
the gradient varies is also important. The designer would 
normally select certain degrees of symmetry. Three are 
most common; axial, radial or spherical gradients. 

GRIN material in the shape of rods now commonly use 
axial and radial gradients in such devices as fiber-optic 
cables and photocopying machines. Spherical gradients are 
not available because spherical gradient material has not 
yet been produced. Spherical gradients do offer an 
advantage in that they lend themselves more easily to 
mathematical description since closed form solutions for 
the ray paths can be obtained. 

A good introduction to gradients and the materials 
currently available for production of large scale GRIN lens 
systems is provided by Moore [Ref. 6], and Light [Ref. 7]. 
These also give the reader an idea of the advantages that 
xmay be realized along with the problems associated with the 
design and production of such items. 

C. AIM AND INTENT OF THIS STUDY 

The aim of this study is to complete the design of the 
Gradient Index Seeker Lens (GISL) initiated by Carr. This 
design uses a spherical gradient and 'pointed' lens shape. 
This study does not determine the aerodynamic efficiency of 
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the design but rather ensures that the lens has a 
generally pointed character that will also allow some 
degree of imaging for off-axis objects. The spherical 
gradient is chosen since it seems to fit the geometrical 
shape of the lens and also since the three-dimensional 
analysis of rays is much simpler. 

This design is accomplished through the use of ray 
tracing techniques. Thus, no actual study of the lens 
material nor the manner in which such a lens could be 
manufactured is intended. It is the hope of the author 
that studies of this nature will provide the impetus for 
the development of suitable GRIN material with spherical 
grad ients . 

Once the lens design was completed a 'best' lens was 
chosen and a simple detector system applied to it to see if 
precise definition of off-axis targets is possible. 
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II. GRIN LENS DESIGN 



A. GRIN LENS DESIGN METHODOLOGY 

The method by which the GRIN Lens was designed in 
Carr's thesis [Ref. 4], is the subject of this section. A 
cross-section view of the top half of the lens is shown in 
Figure 2.1. 

As can be seen, the lens is designed by defining the 
inner surface as a cone and then constructing the outer 
surface so that a backwards ray trace froia a point source 
will result in parallel lines emerging. This is done by 
tracing rays from the back edge of the lens towards the 
front surface and finding the point of intersection of the 
ray and front surface line. The front surface line is 
found by using the previous ray solution as described 
below. The intersection of the ray and the line is 
determined using an iteration procedure. Once a solution 
is found the front surface line that will cause the ray to 
emerge parallel with all other rays is determined. This 
front surface line is then used to find a point solution 
for the next ray. About 1000 rays are traced starting at 
the outside edge of the inner cone and working inwards 
towards the center. Details of the iteration process is 
included in later sections. 
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Once the lens outside surface has been constructed the 



lens is then studied to see how well it will image three 
dimensional skew rays that enter the lens at varying 
degrees of incidence. To accomplish this, a grid is placed 
in front of the lens and then rays are traced from the 
intersection points of the grid through the lens and onto 
the image plane. The grid is then tilted at successive 
angles and the resulting 'spots' on the image plane 
recorded for each case. The manner in which these spots 
move on the image plane and their size are the major 
performance characteristics of the lens. 



B. GRIN RAY PATH DESCRIPTION 



Following Marchand [Ref. 5], GRIN ray paths in 
spherical gradients are found through the solution of the 
appropriate line integral: 



where 




+ e 



(r 



0 



r ( n 



2 2 



r 



dr 

2,1/2 

- e ) ^ 



( 2 . 1 ) 



e = r nr sin = ^ngrgsin q 



( 2 . 2 ) 



and 



' 1 

It should be noted that the value of r is -1.0 if the 
radius is decreasing with angle and +1.0 if the radius is 
increasing with angle. 
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Figure 2.1 GRIN Lens Coordinate System and Parameters (from Carr). 



The choice of index gradient is instrumental in 
allowing the integral to be solved in closed form. In this 
study the index is defined as: 

n2 = (a + br^/R^^ ) 

The constant a defines the strength of the index at the 
origin and the constant b defines the strength of the 
gradient and whether it is negative or positive. 

As can be seen in Figure 2.2, the ray path is described 
by the constant e in conjunction with the parameter i'-, 
the angle between the radial arm of the coordinate system 
and the direction of the ray at any point. 

The solution of the integral in equation 2.1 results in 
the follov/ing expression: 

9 = Oq - j (sin"^(argQ) - sin*"^(arg)^ (2.3) 

where 

a - 2e/rg^ 

z 

and 



arg = 



a - 2e/r 

(a^ + 4be^/R^^)^/^ 



or, expressed in terms of the radius. 
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r 



/2 lei 

{a+ (a^+4be^/R^^ ) sin [-2£(<^9 q) - sin" ^ ( arg^ ) ] 

(2.4) 

C. GRIN LENS GEOMETRY 

A factor which must be considered is the manner in 
which the lens is cut out of the GRIN material. Figure 2.2 
shows an orientation of the material with respect to the 
lens such that the center of symmetry of the material is 
inside the lens. Carr has defined this as a negative value 
of OB, the distance from the point 0 to the point B. 
Carr was able to work out the lens problem for positive 
values of this parameter only, thus the major portion of 
this study is to develop methods that enable solutions for 
negative values of OB. In so doing, two factors have to 
be overcome : 

a) ray matching - When the ray path radius initially 
decreases and then increases with change in angle the ray 
path equation must be matched at the point of minimum 

rad ius . 

b) near-radial lines - Rays that show very little 
curvature are difficult to describe in terras of radius and 
angle and these cause problems in the iterative solution 
scheme . 
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Figure 2.2 GRIN Lens and Ray Geometry. 



The solution to the ray matching problem is relatively 
simple and is discussed in the next section. The 
resolution of near- radial lines is a more complex task and 
is covered in the following chapter. 

D. RAY MATCHING 

As Stated previously the ray path equation changes 
signs as the radius of the path decreases and then 
increases. Such a path is shown in Figure 2.3. Simply 
changing sign in the equation as the point of minimum 
radius is passed does not produce the desired result, but 
rather has the effect shown in Figures 2.3 and 2.4. 

The resolution of this problem is fairly simple. At 
the point of minimum radius the arcsine function in 
equation 2.3 containing this radius (arg) reduces to -/2. 
This is due to the fact that at this point the value of 
is '/2 . Thus: 



sin 4 = 1 



and 



e2 = n^r2 



Thus ; 




a- 2 (a+br^/R ^ ) 
z 



1 



arg 



(a^+4b(a+br^/R^^)r^/R^^) 
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and 

sin” 1 ( arg ) = - tt/2 

If the radius is initially decreasing as shown in 
Figure 2.3, we can find the angle of the minimum radius. 

We call this the critical angle or 9 l: 

0L = Qq “ f {sin'^Carg^) + tt/2 I 

The positive and the negative forms of the equation 
should both contain the point of minimum radius, in other 
words, by approaching the point from the left or the right 
the same result should be obtained. The fact that this 
does not happen is shown in Figure 2.3 and reveals the need 
to match the rays at the point of minimum radius. 

This matching is achieved by simply adding a constant 
to either the positive or the negative form of equation 2.3 
and then matching at the point of minimum radius in the 
follov/ing manner: 

0R = 0L + c 

thus 

- Y { Tr/2+sin”^ (arg^ ) } = + j f Tr/2+sin'’ ^ ( arg^ ) } + c 

The constant c can then be found : 

c = - c I TT/2+sin~^ (arg^ ) } 

This constant must now be applied for all points past 
the point of minimum radius; in our example all points in 
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Figure 2.3 Ray Path Shov/ing Improper Matching and 
Minim.um Radius. 
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Figure 2.4 Ray Path Radius and Angle. 
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the increasing radius portion of the ray must be found 
using the form: 

0 = Sq + j {sin (arg^) = sin (arg)j + c 

In the iteration scheme the ray path equation must be 
more carefully handled. First, the critical angle and the 
constant c must be found. Second, the angle of the 
ray must be checked to see if it has passed the critical 
angle. If it has then the constant c is added to 
equation 2.3 and the point recalculated. The details of 
the iteration methods will be described in the following 
chapter . 

For further clarification Figure 2.5 shows the regions 
in which ray matching is required for a certain position of 
the center of symmetry. Note also that were the value of 
OB positive then all rays in the lens would never reach a 
minimum radius point. Thus ray matching is a problem 
unique to the negative OB case. 
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Figure 2.5 Ray Matching Regions. 
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III. 



NEAR- RADIAL LINES 



A. DESCRIPTION OF DIFFICULTY 

As mentioned previously, the ray paths that are nearly 
radial are difficult to describe mathematically in terms of 
radius and angle. This is due to the simple reason that a 
straight line originating from the center of the coordinate 
system has no change in angle and an undefinable change in 
rad ius . 

The region in which this problem occurs in the lens 
design is shown in Figure 3.1. As can be seen, near-radial 
lines mean that the angle \i; and thus the constant e 
will be very small, so small as to cause equation 2.4 to 
approach a singularity and thus break down. In the 
iterative scheme, as the front surface solutions are found 
from the outside edge of the lens inwards, the values of 
and e become increasingly smaller if the center of 
symmetry is inside the lens (OB negative). When this 
region is reached, equation 2.4 reduces the rg for all 
angles and futher construction of the front surface is 
rendered impossible. This effect is shown in Figure 3.1. 

Of note is the fact that if the radial line region is 
passed successfully, the rays begin to curve in the 
opposite sense, that is they are always convex to the 
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radial lines in the case of a negative gradient and concave 
to the radial lines if a positive gradient is used. 

The resolution of the near-radial lines proved to be a 
difficult task and became a major portion of this study. 

The manner in which it was overcome is the subject of the 
following sections. 

B. PRELIMINARY WORK - IMPROVING THE ITERATION SCHEME 

The first step in the resolution of the problem was to 
simplify the iteration scheme. In his work Carr developed 
a scheme based on the geometry of the rays and coordinate 
system [Ref. 4], which worked well for all positive values 
of OB (that is, for all orientations of the center of 
symmetry to the left of' the inside surface cone. See 
Figure 2.2). It was decided instead to use the Newton- 
Raphson approach, which is simpler in concept and not 
dependent on geometry. 

Carr used the Newton- Raph son method in the second 
portion of the design involving the skew ray trace, but he 
did not use it for the lens construction. 

The manner in which this scheme works is now described 
so that the reader will be able to better understand the 
concept of how all iterative methods work for this design, 
as well as how the ray matching previously described fits 
in . 
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Figure 3 . 1 Ray Trace Showing Region of Radial Lines and 
Effect of Singularity. 
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The geometry associated with this plan is shown in 
Figure 3.2 and a flov/chart outlining the iteration logic is 
shown as Figure 3.3. As can be seen with reference to the 
two figures, an initial point is first chosen by drawing a 
straight line from the inside surface intercept to the 
slope of the front surface obtained from a previous ray 
path solution. Newton- Raph son seeks to find the intercept 
of this front surface line with the curved ray path. This 
is done by using the initial angle to find the radius and 
derivative with respect to angle for each line and then 
calculating a new angle according to the relation: 



e" = 0' - 



r £ - r 

f s gr^ 

• • 

r^ - r 
f s gr 



This new angle is then used to recalculate the radii 
and derivatives and thus an iterative loop if formed. Once 
the difference in angle is small enough tlie iteration is 
stopped and the solution considered to be found . 

Ray matching is also employed in this scheme. As each 
new angle is calculated, it is checked to see if the 
critical angle has been passed. If it has, the constant c 
of equation 2.5 is added to the angle and thus the proper 
radius and derivative will be then found. 

Once a solution has been found it is used to find a new 
front surface slope by determining that slope which will 
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Figure 3.2 Newton-Raphson Iteration Geometry. 
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Figure 3.3 Newtor>- Raphson Logic Flowchart. 
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cause the refracted ray to emerge parallel to the axis. It 
is this slope that is used to find the solution for the 
next ray in the scheme. 

When this new iterative procedure was implemented it 
was found to work efficiently for both the positive and 
negative gradients, but it did not resolve the problem of 
near- radial lines. 

The reader should refer to Appendix A for a listing of 
program GISL, which is the ray tracing program originally 
written by Carr and modified by the author. 

C. PRELIMINARY WORK - DOUBLE PRECISION 

The next avenue that was explored was to place the ray 
tracing scheme in double precision. This was considered 
necessary since very small differences in angle were 
required in the region of the radial lines. 

Once this was installed in the routine it was found to 
improve significantly the accuracy of the calculations and 
succeeded in delaying the point at which the iterations 
broke down. It did not solve the crux of the difficulty 
however . 

D. PRELIMINARY WORK - CHECKING ACCURACY 

At this point it was decided to develop a scheme to 
check the accuracy of the front surface solution. The 
manner in which this was done is shown in Figure 3 .4 
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The ray trace was completed in a forward direction froia 
the front surface to the focal point in a two-dimensional 
fashion. These rays were caused to originate from the mid 
point of each set of points that define the front surface. 
The slope of the front surface for each of these rays was 
taken to be the line drawn between the two adjacent points 
on the front surface. In this fashion the worst case v;as 
chosen since the skew ray trace used to determine the 
performance of the lens uses incident rays that can hit the 
front surface at any point. 

These rays were then traced through the lens and back 
towards the focal point. The x- coordinate of the focal 
point was defined and the y-coordinate calculated from the 
ray trace to see how close it would come to the 
y-coordinate of the focal point originally used to start 
construction of the front surface. A listing of this 
routine is shown as Appendix B. 

Another benefit of this two-dimensional ray trace was 
that it allowed the author to confirm the methods by which 
a forwards ray trace is done. In this manner ray matching 
and sign convention in angles was checked. Also this ray 
trace is different in that a solution between the ray path 
in the lens and the back surface, rather than the front 
surface, was now required in the iterative solution scheme 
and this could now be verified to ensure that it was 
properly done. 
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Another method of checking the accuracy of the solution 
would be to study the front surface smoothness in the 
region of near- radial lines. Although not as precise as 
the above method this idea allowed the author to confirm 
that the front surface was sufficiently smooth. By this 
means it was also possible to see a physical representation 
of the lens in the region of the singularity. 

E. FIRST PROPOSED SOLUTION - STRAIGHT LINES 

The first proposed method of solving the mathematical 
difficulty was to simply replace the near- radial lines with 
completely straight lines just before the singularity was 
reached . This could be done by following the decrease of 
the parameter i|;^ until it became small enough to ensure 
that the lines were nearly straight. 

This idea was implemented but it was found not to be 
sufficiently accurate. It must be emphasized that an 
inaccuracy in even a single point could not be tolerated 
since this would mean that all succeeding points in the 
iterative solution scheme would carry that error as well. 

For these reasons the straight line solution was 
rejected and other means of solving the difficulty were 
sought . 

F. SECOND PROPOSED SOLUTION - REDUCED EQUATIONS 

The second method that was used to resolve the problem 
of near-radial lines was to reduce equation 2.4 by 
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start of trace 
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employing small angle aproximations , binomial expansions 
and curve fitting. The details of how this was done 
follows . 

Recalling equation 2 . 3 : 

e = Gq - f {sin'^Carg^) - sin'^(arg)^ 

The arguments of these arcsine functions were reworked 
as follows; 



and 



where 



argo 



arg 



a-2 (a+b ' ) <1^0^ 

( a^+ 4 b ' ( a+b ' ) 'I-q ^ 

a-2 (a+b ' ) ii-^^r^^/r^ 
(a^+4b'(a+b')i)';o^)^'^^ 



1^0 = sini[;Q 

b' = br02/R22 

Note the angle i;;0 is very small. 

These arguments were then reduced using a binomial 
expans ion ; 

2(a+b’)^fe^ 
argo = ^ 2 
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Figure 3.5 Expanded View of Front Surface - Straight Line 
Method . 
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and 



arg 



1 - 



If 

2(a+b')'!'o^(-§- a+b') 



Again the i|)q term ensures that the expansion is valid. 

A function of the type arcsin(x) where x is known to 
be close to the value 1 can be curve fitted into the form: 

sin"l(x) = tt/ 2 - [2(1 - x)]l/2 



This equation applies to both arcsine functions as the 



0 / 



0 



term ensures that the argument is close to 1. 

When placed in the equation the arguments result in the 
following final expression: 

2 



f. 



e - 



° ((a+b’) - [(a+b')( — 



a+b’ 1 



In terms of radius: 



^ ^ ^ (a(a+b* ) 

([-^ (0 - 6.) - (a+b')J^ - b' (a+b- ) 

0 

The reduced equation was then implemented in the 
iterative front surface solution scheme by constantly 
checking it against the complete equation. When the 
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difference between the two became small, the reduced 
equation was subsequently used in solving for the front 
surface. The reduced equation was then turned off when the 
rays had progressed an equal distance on the other side of 
the singularity. See Figure 3.2 for clarification. 

The equation worked well in resolving the front surface 
for the lens construction and did produce a smooth surface 
similar to that shown in Figure 3 .5 . 

It did not do well when the accuracy program check was 
applied. One reason is that although it was easy to 
implement the equation for the lens construction it was not 
so easy to complete the forwards trace since these rays hit 
the front surface randomly making it difficult to know when 
to use the reduced equation. In practice the reduced 
equation could be implemented when the absolute value of 
the parameter \};q became smaller than that used in the 

lens construction phase during the time that the reduced 
equation was used. There was then no way of telling how 
accurate the front surface was, in other words, there was 
no method to show how the rays would physically react in 
the region of the singularity. Thus the reduced equation 
method would work but it is not an entirely satisfactory 
solution . 
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Figure 3.6 Reduced Equation Impleinentetion . 



Program CHECK was tried without using the reduced 
equation but disastrous results were obtained for those 
points in the region of the singularity. 

For these reasons it was decided that another means of 
resolving the problem was required. 



G. THIRD PROPOSED SOLUTION - MARCHAND METHOD 

The third proposed solution is drawn from Marchand 
[Ref. 5J p. 27, in which he proposes a means of improving 
the ray trace formula so that nearly radial lines will be 
more easily described . 

The method is based on redefining variables as 
follows ; 

0 = eMo (3.1) 

where 



«0 = J 



r dr 

>^0 r(n2r2 - 



Marchand describes the ray path using the parameters 
and B such that; 



3 = sin9/sin>j;g 



(3.2) 



a = cose - BCOSibQ 

This enables the coordinates of the ray to be described 



as : 



X 



r(a— + 3—) 
r* n 

0 0 
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Y 



= r ( ( 



^0 



- 20 , 

^0 



z 








Upon substitution of equations 2.2 and 3.2/ equation 
3 .1 becomes : 



^ = rnQrQMo sine" 

The parameter can now be more accurately solved for 

small angles of ') using a Maclaurin series for the sine 
function . 

This method was then implemented in program GISL using 
an iteration scheme based on Cartesian rather than polar 
coord ina tes . 

This method was not found to improve the accuracy of the 
ray trace in the region of near radial lines and was not 
adopted . 

One reason for this may lie in the fact that in our 
solution to the integral in equation 2.1 the value e no 
longer appears at the front of the expression. Thus there 
is no advantage to redefining the variable as shown in 
equation 3.1. 



H. FOURTH PROPOSED SOLUTION - ANGLE IN TERMS OF RADIUS 
The fourth method and the one that proved successful 
involved the solution of the front surface slope and the 
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GRIN ray using the form of the ray trace equation shown in 
equation 2.3 rather than that of equation 2.4. 

The reason for this choice is shown in Figure 3 .7 . 

This figure shows that when the rays are nearly radial a 
small change in angle will result in a huge change in 
radius. Thus the form of the ray tace equation using r( ) 
is inherently unstable. Previously we have used this form 
in order to make use of the Newton- Raph son technique. It 
was thought that if the iteration scheme could be reworked 
so that the form 0(r) of the equation could be used then 
better results might be obtained. 

A suitable iteration method was found and is 
illustrated in Figures 3 .8 and 3 .9 . 

The figures show that an initial point was chosen in 
the same manner as before, that is by drawing a straight 
line from the ray refracted at the inner surface out to the 
front surface slope. This point provides the radius for 
equation 2 .4 so that a new angle can be calculated . This 
angle is then used to find the corresponding radius out to 
the front surface slope. This new radius is then used in 
equation 2.4 to find the next angle and thus an iterative 
loop is formed . 

Thus a simple scheme was found . It should be noted 
however that this method will only converge when the 
initial angle is greater than 90 degrees. The reason for 
this lies in the geometry and is shown in Figure 3 .10 . 
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Figure 3.7 Radius and Angle Relation. 



Thus the iteration will only converge v;hen the rays are in 
the region of the nearly radial lines. This does not cause 
a difficulty since the Newtor>- Raphson method, modified to 
include ray matching, can find solutions efficiently in the 
region where the initial angle is less than 90 degrees. 

Figure 3.10 is also drawn to show the reader an example 
of a forward ray trace. 

This simple scheme was implemented and found to work 
successfully for both the lens construction and the skew 
ray trace; see [Ref. 8]. 
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Figure 3.8 Angle in Terms of Radius Iteration Logic 



Flowchart . 
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Figure 3.9 Angle in Terms of Radius Iteration Geometry. 
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Figure 3.10 Diverging Iteration - Forward Ray Trace. 
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IV. GRADIENT INDEX SEEKER LENS RAY TRACE RESULTS 



A. FRONT SURFACE ACCURACY 

As mentioned previously, the front surface of the 
Gradient Index Seeker Lens was checKed for accuracy by 
completing a ray trace from the midpoint of each set of 
solutions. A sample of the output from program Check is 
shown in Table I so that the reader will have an indication 
of the accuracy of the ray trace in the region of near- 
radial lines. Note that the most radial lines occur when 
the absolute value of the angle \||q is smallest. 

B. SPOT SIZE RESULTS 

The results of program GISL showing the variation of 
spot size are shown in Figure 4.1. This figure shows how 
the spot size varies as gradient strengths and orientations 
of the center of symmetry are changed for a median angle of 
incidence of 0.3 radians. 

Note that all orientations of the center of symmetry 
are now included. The studious reader should refer to 
Carr's work [Ref. 4], to compare results. 

C. SELECTION OF A BEST LENS 

The best lens was selected by determining those 
parameters that would result in a small spot size and also 
cause the least number of rays to be reflected or refracted 
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outside the missile diameter. Figure 4.1 shows that these 
criterion were not easily discernable and a judgemental 
decision was required. 

The reader should note that a value of 2.25 was chosen 
for the parameter a. This was done since the resulting 
value of 1.5 for the index of refraction at the center of 
symmetry is common to a large number of possible lens 
materials. Also of note is the fact that the focal point 
was chosen to be 2.0 units along the x-axis from the inner 
surface cone. This was done since this value resulted in 
smaller spot sizes generally. 

With these criterion in mind a value of -0 .2 for 08 
and a 5% negative gradient were chosen as parameters for 
the best lens design. 

A more detailed study of the lens parameters could be 
accomplished with use of an optimization program in which a 
large number of parameters could be varied in order to find 
the best lens. This would achieve a more satisfactory 
result than the one shown. 
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TABLE I 



FRONT SURFACE ACCURACY 



1 Line Number 


1 Angle ’.Vq 
1 (Radians) 


1 Error in y- I 
1 coordinate I 


number of 
iterations 




1 895 


~i 

1 0.0411 


1 r 

1 0 .0008 1 


4 




1 896 


1 0 .0390 


1 0 .0007 1 


4 




1 897 


1 0 .0369 


1 0 .0008 1 


4 




1 89 8 


1 0 .0349 


1 0 .0009 1 


4 




1 899 


1 0 .0328 


1 0 .0007 1 


4 




1 900 


1 0 .0308 


1 0 .0009 1 


4 




1 901 


1 0 .0288 


1 0 .0008 1 


4 




1 902 


1 0 .0268 


1 0 .0008 1 


4 




1 903 


1 0 .0247 


1 0 .0009 1 


4 




1 904 


1 0 .0227 


1 0 .0009 1 


3 




1 905 


1 0 .0207 


1 0 .0008 1 


3 




1 906 


1 0.0137 


1 0 .0009 1 


3 




1 907 


1 0.0167 


1 0.0008 1 


3 




1 90 8 


1 0.0147 


1 0 .0008 1 


3 




1 909 


1 0.0127 


1 0.0010 1 


3 




1 910 


1 0.0107 


1 0 .0009 1 


3 




1 911 


1 0 .0087 


1 0 .0008 1 


3 




1 912 


1 0 .0067 


1 0.0010 1 


3 




1 913 


1 0 .0048 


1 0 .0009 1 


3 




1 914 


1 0 .0028 


1 0.0010 1 


2 




1 915 


1 0 .000 8 


1 0.0010 1 


2 




1 916 


1 - 0.0010 


1 0.0010 1 


2 




1 917 


1 -0 .0030 


1 0.0010 1 


3 




1 918 


1 -0 .0050 


1 0.0010 1 


3 




1 919 


1 -0 .0069 


1 0 .0009 1 


3 




1 920 


1 -0 .0089 


1 0.0011 1 


3 




1 921 


1 - 0.0108 


1 0.0010 1 


3 




1 922 


1 - 0.0127 


1 0.0010 1 


3 




1 923 


1 - 0.0147 


1 0 .0009 1 


3 




1 924 


1 - 0.0166 


1 0.0011 1 


3 




1 925 


1 - 0.0185 


1 0.0010 1 


3 




1 926 


1 -0 .0204 


1 0.0011 1 


4 




1 927 


1 -0 .0223 


1 0.0011 1 


4 




1 928 


1 -0 .0243 


1 0.0010 1 


4 




1 929 


1 -0 .0261 


1 0.0012 1 


4 




1 930 


1 -0 .0281 


1 0.0011 1 


4 




1 931 


1 -0 .0300 


1 0.0012 1 


4 




1 932 


1 - 0.0318 


1 0.0012 1 


4 




1 933 


1 -0 .0337 


1 0.0012 1 


4 




1 934 


1 -0 .0356 


1 0.0011 1 


4 




1 935 


1 -0 .0375 

1 


1 0.0013 1 

1 1 


4 
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V. MIRROR AND DETECTOR SYSTEM 



A. DESCRIPTION OF SYSTEM 

A scanning mirror and a detector were added to the GRIN 
best lens in order to determine if precise definition of 
off-axis targets would be possible. The design of these 
two new elements is shown in Figure 5.1. 

The mirror is made to 'nod' at a high rate and thus 
scans the spot across the detector. The detector consists 
of two sections, as shown in Figure 5.1. Thus, when the 
spot is scanned across the two plates a signal from each 
detector will be generated and will result in the form 
shown in Figure 5.2. This arrangement will allow an 
accurate measurement of the mirror angle for which the spot 
crosses the axis. The electronics of the detector signal 
processor could then be arranged to relate the mirror angle 
to the angle of incident radiation so that the precise 
1 ine-of- sight angle of an off-axis target will be known. 

Note that this is a simple detector system that will 
only resolve angles in one plane. A more sophisticated 
system along the same lines could be designed with a four- 
element detector and a more intricate method of scanning 
the spot across the detector. It was decided to keep the 
mirror and detector system as simple as possible since the 
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aim of this work is only to show that a precise angle 
measuring device is feasible. 



B. LENS TO DETECTOR RAY TRACE 

This section describes the mathematics of the ray trace 
from the lens to the mirror and then to the detector. A 
listing of program DETPLOT which accomplishes this task is 
included as Appendix C. 

The ray trace consists of the following parts; 

a. ray direction exiting lens 

b. point of intersection on mirror 

c. ray direction exiting mirror 

d . detector ray count 

1 . Ray Direction Exiting Lens 

The ray direction exiting the lens is derived from 
the direction cosines of the ray as calculated by program 
GISL. This data is used as an input to program DETPLOT 
along with the coordinates of the point at which the ray 
left the lens. 

2 . Point of Intersection on Mirror 

The point of intersection of the ray and mirror may 
be described by the following expressions; 



Xm = DojnK + XQ 



(5 . 1 ) 



Ym = D.niL + yg 



(5 . 2 ) 



z 



m 




( 5 . 3 ) 
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Figure 5.1 Lens, Mirror and Detector System. 
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Figure 5.2 Detector Signal. 
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MIRROR ANGL 



since the surface of the mirror is tilted along the 
z-axis it can be described by the equations: 

= tan9^vi(Xf^ - MP) (5.4) 

Yi^2 cos2 9 ^ + Zj^2 ^ r (5 .5 ) 

Where equation 5 .4 describes a tilted plane and 
equation 5.5 limits the plane to the mirror diameter. 

Neglecting equation 5 .5 for a moment, equations 
5.1, 5.2 and 5.4 may be combined to solve for D^j^. The 
result is: 

^ ^ tan e^(xp - MP) - yg 

^?.m L “ K tan 6 

m 

Now that is known, equations 5.1, 5.2 and 5.3 

may be used to find the coordinates of the mirror 
intersection point. These are then checked with equation 
5.5 to make sure they are not outside the mirror diameter. 

3 . Ray Direction Exiting Mirror 

The solution of the ray direction leaving the 
mirror is best understood vectorially. See Figure 5.3. 

Reflection off a mirror is governed by two factors. 
Firstly, the angles of incidence and reflection must equal. 
Using vectors this is expressed as: 



->• 

n 



r . 
1 



= - n 



r 



In component form 



nr . + n r . 
X XI y yi 



n r 
X xr 



n r 

y yi: 



(5 .6 ) 
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Secondly, the ray is reflected in the same plane as 
the incident ray. This may be expressed vectorially by 
using the cross product to generate the vector np as 
shown in Figure 5.3. Thus: 

rip = r^xn = r^xn 
In component form: 



ex(-nyr 2 i) + ey(-nxrzi) + ez ( nyr^i- n^ry ^ ) = 



-y -y 

®x(”*^v^zr) ®v(“^x^zr) ■*" ®z(*^y^xr “ *^x^vr) 



'X \ “y'-zr^ ^ '^y V "x*-zr 
Three equations result: 



‘x^yr 



" *^y^zi 


- -nyrzr 






= "*^x^zr 




*^x^y i 


= nyrxr - nxry^ 


(5 .7 ) 



"y‘-xi 

The first two results are not unexpected as the 
mirror should not effect the z- component of the ray. 

Equations 5 .6 and 5 .7 may now be combined to find 
the X- and y-components of the ray. When this is done the 
two following expressions result: 



/ 2 2 , 
(n - n ) 
y X 



2n n 

X y 



xr 



■yr 



(n^ 

' X 


+ 


2 

n 

y 




(n^ . 

X 


2 , 
"y > 


/ 2 
' X 


— 


2 

n 

y 


) 

y* — 






/ 2 
' X 


+ 


2 

n 

y 




(n2 + 

^ X 


) 

y 



yi 



XI 
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Figure 5.3 Reflection of Ray on Mirror. 
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Since the vector n is entirely in the x-y plane, 
its components may be related to mirror angle as follows; 

nx = -sin^ni 

ny = cosBni 

When the above equations, along with well known 
trigonometric identities are used the following expressions 
will be found . 

rxr = cos(2%)rxi + sin(2%)ryi 

ryr = -cos(2 ‘3m)ryi + sin (2 '^ 1 ^) rxi 
Thus these equations represent the change in direction for 
the exiting ray and form two of its direction cosines, the 
third being 

4 . Ray Count on Detector 

The coordinates of the ray at the detector are 
described by the expressions: 

Xd = Dnd K' + Xm 

Yd = Dnd L' + Ym 
zd = Dnd + zm 

Since x^^ is used to define the position of the 
detector, it is known. Thus the value may be solved 

and the corresponding y and z-components may be found. 

The ray is then checked to make sure it is within 
the detector diameter. If it is the y-coordinate is 
checked to find out if it is negative or positive, and thus 
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whether the ray has struck the upper or lower half of the 
detector. The ray count is then incremented accordingly. 

Program DETPLOT traces all rays for each increment 
of mirror angle and thus generates plots such as that shown 
in Figure 5.2. 

C. MIRROR AND DETECTOR SYSTEM RESULTS 

Appendix D shows the results of the detector signals as 
a function of mirror angle at angles of incidence from 0 to 
40 degrees. For these the mirror pivot point was placed at 
1.4 units along the x-axis from the inner surface cone 
since this was the smallest distance that would still allow 
full movement of the scanning mirror. The focal length was 
adjusted to 2.8 units for the lens design calculations so 
that the spot would be imaged as small as possible on the 
detector. Figure 5.4 summarizes the results shown in 
Appendix D. As can be seen, the smaller angles of 
incidence (those in the range 0 to 16 degrees) result in a 
linear relation between angle of incidence and mirror tilt 
angle, corresponding to an on-axis spot. The larger 
values of the angle of incidence resulted in a more random 
relation between these two quantities. 

Appendix D also shows that the error in determining 
mirror angle increases as the angle of incidence increases. 
This is due to the fact that each successive plot is less 
distinct and thus a precise measure of the angle becomes 
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more difficult. Although this effect shows a deterioration 
in the lens system performance, it is not harmful since, as 
the missile corrects itself, the measurement of angle will 
become more accurate. 

For these reasons this simple system is considered to 
be a viable means of target detection and shows that the 
GISL can be a successful missile seeker lens. 
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Figure 5.4 Mirror Tilt for On- Axis Spot and Angle of 
Inc idence . 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

Two objectives were accomplished in the course of this 
study. Firstly, the Gradient Index Seeker Lens design, as 
initiated by Carr [Ref. 4], was completed in that all 
orientations of the center of symmetry with respect to the 
lens were made possible to describe mathematically. The 
results obtained for the spot size on the image plane 
showed trends that indicated the parameters that would 
produce a 'best' lens, but precise definition of these 
values was not possible. 

Secondly, it was shown that a scanning mirror and 
detector system applied to the lens could be described 
using ray- trace techniques and also that such a system 
could be used to provide precise angle measurements of 
off-axis targets at angles less than 16 degrees. For 
larger values of angle, error in definition became 
increasingly greater. 

Specifically, then, the following conclusions are made 

a) The Gradient Index Seeker Lens can be described 
using ray- trace techniques and can be shov/n to provide 
adequate imaging of off-axis objects. 

b) A scanning mirror and detector system can be used 
to provide precise angle measurements. 



68 



B. RECOMMENDATIONS 



The following recommendations for further study are 
made : 

a) An optimization program should be applied to 
program GISL to study the variations of all important 
parameters . 

b) A detector system capable of detection in all 
directions should be designed. 

c) The manner in which signal processing of the 
detector signals can be carried out should be studied. 

d) A study of Gradient Index material and the methods 
by which a lens could be manufactured should be undertaken. 
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Figure I>4 . Detector Signal at 16 Degrees Incidence. 
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Figure D-6. Detector Signal at 24 Degrees Incidence. 
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